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SUMMARY 

Supersonic flow against blunt bodies placed in boundary layers or 
wakes is discussed qualitatively, and it is concluded that wedge-shaped 
or conical dead-air regions should form ahead of the body if part of 
the upstream velocity profile is subsonic and if the body is fairly 
thick relative to the initial boundary layer or wake thickness. Quan- 
titative analysis of this type of flow indicates, however, that for 
each free- stream Mach number, there is also a maximum relative body 
thickness beyond which wedge-type or conical-type separation regions 
cannot occur. This maximum relative body thickness is large for high 
Mach numbers but approaches zero for a Mach number of 1. 

For the intermediate range of body thicknesses, an analysis of the 
two-dimensional flow against blunt bodies mounted on a flat plate 
agreed quantitatively with experimental results in the Mach number range 
of 1.73 to 2.02. As the body thickness approached the maximum theoreti- 
cally possible for wedge-type separation, unsteady flow was observed. 

For bodies of the same order of thickness as the initial boundary 
layer, experimental pressures near the leading edge of the bodies 
remained higher than those predicted for wedge-type separation. 


INTRODUCTION 

When a blunt body is placed in a uniform supersonic stream, a 
detached shock wave forms ahead of it. This detached shock is normal 
in front of the body and decays gradually into a Mach wave at large 
distances from the body. If, however, the blunt body is placed in a 
region of nearly constant static pressure, but nonuniform stagnation 
pressure, such as a boundary layer or a wake, it becomes impossible to 
satisfy certain flow requirements unless the form of the shock is 
altered. Consider, for example, the sketch in figure l(a) of a blunt 
body in a completely supersonic wake. If the static pressure ahead of 
the shock is assumed to be constant and the velocity gradient to be as 
indicated, the pressure behind a normal shock at the stagnation stream- 
line would be less than the pressure behind the shock in the free stream. 
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With the dashed shock form shown, there would consequently be a posi- 
tive pressure gradient toward the center line and a convergent flow 
behind the shock, which is an impossible situation. In any real flow, 
the pressure gradient normal to the stagnation streamline on the axis 
of symmetry must be zero or negative. Since the Mach number at the axis 
is lower than in the free stream, no shock will produce pressures at the 
axis as high as those behind the shock outside the shear layer as long 
as the stagnation point remains on the body. The required high pres- 
sure near the vertex of the shock can be attained only if the stagna- 
tion point moves upstream^ hence a dead-air region will form (fig. 1(a)) 
which effectively changes the nose of the body to a form compatible 
with the flow requirements. 

When the velocity near the axis is subsonic, as in a boundary 
layer, the high pressures in the dead-air region can feed upstream and 
modify the flow outside the boundary layer or wake until some sort of 
equilibrium is established. If the body is fairly thick relative to 
the shear layer, the dead-air region can become large relative to the 
thickness of the shear region. If the circulatory motions in this 
region are negligible, the boundary of the dead-air region must be a 
const ant -pres sure surface. Hence, the form of the dead-air region out- 
side the boundary layer or wake should be that of a wedge for two- 
dimensional flow and a cone for axially symmetric flow (fig. 1(b)). 

As the body thickness is increased indefinitely, it seems reason- 
able that a configuration different from that shown in figure l(b) will 
occur. When the shear layer becomes insignificantly thin relative to 
the body, the detached-shock configuration encountered in uniform 
supersonic streams should reappear. 

The purpose of the present paper is to investigate the nature of 
the equilibrium that determines the size of a wedge- type or cone-type 
separation region and the limits beyond which this type of separation 
cannot exist. For two-dimensional flow against blunt bodies mounted 
on a flat plate, theoretical predictions are compared with experimental 
results. This investigation was conducted at the NACA Lewis laboratory. 


ANALYSIS 

The initial analysis will be concerned with the flow against two- 
dimensional bodies mounted on a flat plate. The symbols used are 
defined in appendix A. 

The analysis is based on the simplified picture of the separation 
phenomenon shown in figure 2. The curvature of the shock resulting from 
boundary- layer - shock-wave interaction is neglected, and the shock is 
assumed to remain straight and of constant intensity from the edge of 
the boundary layer to the point where the line perpendicular to the 
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dead-air "boundary from the corner of the body intersects the shock. 

The dead-air region is assumed "bounded by a straight line inclined at 
the angle X to the flat plate and to have a constant pressure p^ 

equal to the pressure behind the oblique shock. These assumptions 
should be approximately valid if the initial boundary-layer thickness 
5 q is small compared with the obstacle thickness b and if the circu 
latory motions in the dead-air region are negligibly small. 


The continuity equation for the configuration in figure 2 is 


P&r 


p o u oW + 


pu dy = P-^ 7 !" 6 !) + 


pu ds 


which can be converted to 


Vo (y oV> - c i u i (y r 5 i‘ > 


(i) 


Since there is no friction force along the dead-air region, the 
momentum equation reads: 

(p i" p o )y o + P 1 U 1 2 cos X W ' p oV ( W + 


P5, 


0 


pu' 


^ cos \ ds - 


0 


pu^dy = 0 


0 


which, when combined with equation (l), can be converted to 

,2 


n u i x p i" p o 

y n 1 - — cos X p 

0 1 o P 0 V 


9 0 ‘ 


Pl u l 

p 0 




( 2 ) 


where © and 9-^ are the momentum thicknesses ahead of and at the 
end of the dead-air region, respectively. From oblique shock relations, 


U -, 
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p l" p 0 1 
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so that equation (2) becomes 


where 


and 


©1 1+2 


e. 


Q. 
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Pl u l 


cos 


X 




V 

P 0 Mq 2 


COS 



The parameters P.^, <p, and X are related by the equation 


tan X = 


\ Pi cotep 



( 3 ) 


Hence equation (3) determines the equilibrium values of shock angle, 
detachment angle, or pressure in the dead-air region if 9]_/9 q I s 

known. Precise determination of 0p depends on a knowledge of the 
profile between the dead-air region and the outer stream in the region 
downstream of the shock wave. Since this knowledge is not presently 
available, it will be of interest to see if a valid hypothesis concern- 
ing the magnitude of ®i/9q can ma( ^ e " 

The sketch of figure 2 suggests that a good estimation of ©p 
might be obtained by assuming that the increase of © between 
station 0 and 1 due to mixing is proportional to the increase that 
would be expected if a solid wedge with angle X replaced the dead-air 
region. With this assumption, the value of 9p becomes 


© 
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n l 


9 o + 


e T i 

2 

Pl U l 


dx 

cos X 


where c is a factor of proportionaltiy. 
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When an average shear stress is used for the wedge, this equation 
becomes 


9 


1 


£T l,av 

2 

Pl u l 


_b 
sin \ 


(4) 


With the assumption that the average shear stress on the replacing 
wedge differs little from the shear stress on the flat plate immediately 
ahead of the wedge, equation (4) becomes 


9 1 " 9 0 


Pl U l 


0 sin \ 


because 


-I°_ m M m9 , 

2 l dx / , W 0 

Vo v /x - i 


Rearrangement of the terms yields 


^1 P 0 U 0 €0 o' b 

0 n “ 1 + _ 2 0 sin \ 

0 o 


(5) 


Equating the values of 0 ^/©q from equations (3) and (5) yields 


or 


P 0 U 0 2 eG o' b 1 + 2 h 0 P l 

1 + 2 © sin A “ „ 2 

PpU! 0 Pp 11 ! 


cos \ 


p o u o 


0 * 

£ -- b a (1 + luP-. - Q) tan \ 3 T tan \ 
y 0 ^ 


( 6 ) 


Wow, both for turbulent and for laminar flows along flat plates, 
formulas exist for the skin friction in the form 
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- k(Mo) 
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(7) 


where T] is 2 for laminar flow and approximately 7 for turbulent flow. 
Integration of equation (7) yields 

e o - e ' 1 V 1/n 

Jo 


Hence , 


ft ' 

y o 


_ HiL 

T]Z 


( 8 ) 


and is independent of Mach number and Reynolds number to the extent 
that T] is independent of these parameters. 

Equation (6) now becomes 

€ ^ = — L T tan \ ( 9 ) 

l T)-l 


b m T tan \ (io) 



since 


l s L - b ctn X 


When e is specif ied, equation (10) gives the theoretical vari- 
ation of 9, \, or P-^ as a function of the thickness of the obstacle 
if fHe boundary layer ahead of the separated region is either completely 
turbulent or completely laminar. It appears reasonable that, for 
laminar boundary layer, the effective shear over the wedge should 
nevertheless be assumed equal to that of a turbulent boundary layer, 
inasmuch as the mixing region is probably turbulent regardless of the 
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state of the initial boundary layer. This modification is discussed in 
appendix B, where it is shown that the variation of b/L with \ can 
be calculated quite simply for this case from the values obtained for 
purely laminar flow. Modifications required when the boundary layer is 
partly laminar and partly turbulent are discussed in appendix C. 


For turbulent boundary layers, values of h^ required in equa- 
tion (10) have been calculated as a function of Mq in reference 1 for 

U.-/ zV /N 

“0 l 8 J ‘ 


various power profiles of the form 


For laminar boundary 


layers, h Q can be obtained from the profiles derived in reference 2. 
For the case of zero temperature gradient and a Prandtl number of 0.72 
the equations of reference 2 yield 


1.73 + l.ll(r-l)M 2 
0.664 


(H) 


Values of h as a function of M, obtained from reference 1 for 
N = 7 and from equation (ll) are plotted in figure 3. With these 
values, the variation of X and 9 with b/L for turbulent and laminar 
boundary layers was computed for several Mach numbers for € » 1.0. The 
results are plotted in figure 4. 1 It will be noticed from this figure 
(particularly from fig. 4(b)) that, for each value of Mq, a maximum 

value of b/L is obtained beyond which the present analysis yields no 
solution. This maximum value occurs when the boundary-layer separation 
angle reaches the maximum flow-deflection angle for an attached oblique 
shock wave. For Mach numbers close to 1, wedge-type separation appears 
to be possible only for small values of b/L. For such small values of 
b/L, however, the effects of the boundary- layer shock-wave interaction 
become important. Hence, the present analysis is probably of interest 
primarily for high Mach numbers, where wedge-type separation is pos- 
sible for fairly large values of b/L. 

For values of b/L greater than the maximum value possible for 
wedge-type separation, it was expected that one of two possible con- 
figurations might occur: (l) A configuration sim i lar to figure l(a) 

might form ahead of the obstacle; or (2) the separation wedge might 


"TThe effect of using other values of € , as can be seen from equation (lO)j 
would be to shift the curves of figure 4 toward larger values of b/L for 
€ < 1 and to smaller values of b/L for € > 1 . No estimate of the mag- 
nitude of e can yet be made from theoretical considerations, although it 
appears likely that the momentum change along a free surface should be 
less than along a solid surface. However, details of the shock - boundary- 
layer interaction may invalidate this expectation. 


8 


MCA IN 2418 


extend to the leading edge of the plate. Actually, as will be shown 
later, neither of these configurations was observed. Instead, the flow 
became unsteady even before the maximum analytical value of b/L was 
reached. 

For large values of M^, the present analysis yields negative 

values of b/L for weak shock waves. These negative values probably 
indicate again that details of the shock - boundary- layer interaction 
cannot be neglected when b is of the order of the thickness of the 
boundary layer. It will be noticed from equation (4) that the values 
of cp or X obtained at b = 0 are those which would result if 
were assumed equal to 9 q. With such an assumption there would be no 
variation of separation angle with body thickness, and separation 
angles other than zero would be obtained only at high Mach numbers. 


APPARATUS AND PROCEDURE 

A preliminary experiment to check analytical results was conducted 
in the NACA Lewis 18- by 18-inch tunnel, wherein a series of bars of 
rectangular cross section were mounted on a flat plate with a span of 
16 inches. The angle of attack of the flat plate was variable such 
that a r ang e of free-stream Mach numbers from 1.735 to 2.03 was obtain- 
able. The Reynolds number per foot at a Mach number of 1.9 was 
3.24X10 6 . 

For most tests, transition of the boundary layer was artificially 
induced by placing a l/4-inch strip of carborundum dust across the 
plate, 3/4-inch from the leading edge. This strip was removed for some 
of the tests to check the effect of changing the initial velocity 
profile • 

The flat-faced rods were mounted with their leading edges 8^ inches 

from the leading edge of the plate. For values of b up to 3/4 inch, 
the rods extended over the full span of the plate. For thicker rods, 
however, choking of the tunnel occurred with full-span rods. Some of 
these rods were therefore reduced in length to obtain data for the 

larger values of b/L. The l— inch thick rod, which was the thickest 

for which steady flow was obtained, had to be reduced to a span of 
6 inches to avoid choking the tunnel. The effect of such reduction in 
span is shown in figure 5 for b s 0.5 inch. If it is assumed that the 
extent to which the flow is two-dimensional depends chiefly on the 
span- thickness ratio of the body, the results plotted in figure 5 

indicate that for li-inch rod ( span- thickness ratio » 4.8) the boundary- 
layer detachment angle was reduced about 1*^ below that which would have 
been obtained for a more nearly two-dimensional body. 
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The variation of 0 Q with l along the flat plate with artificial 

transition was obtained at Mach number 1.88 in a previous investigation 
and is shown in figure 6(a) . For comparison, a curve represented by the 
empirical formula 


®0 - 1 “(“o' 1 V 1/7 

is also shown where k(M^ ) was taken equal to 0.0104. This value of 
k(MQ) was obtained from the equations of reference 1 for a Mach number 

of 1.88. The empirical curve is somewhat higher than the experimental 
values probably because transition was induced artificially 3/4 inch 
from the leading edge. The boundary- layer thickness, according to 
reference 1, is about 12 times the momentum thickness at a Mach number 
of 1.9. In figure 6(b), values of 0 o '/ e o obtained from the experi- 
mental data and from equation (8) are plotted. The good agreement 
indicates that little error in the variation of b/L with separation 
angle should result from the use of empirical friction formulas. No 
boundary -layer surveys were made with natural transition on the plate. 


COMPARISON OF THEORY AND EXPERIMENT 

Schlieren photographs of the flow past the rectangular rods are 
shown in figure 7. The waves arising from the transition strip (which 
can be seen in most of the photographs) did not affect the determi- 
nation of the Mach number upstream of the separation wedges since this 
Mach number was calculated using the static pressure measured by an 
orifice located l/2-inch downstream of the transition strip. Another 
weak wave, about two-thirds of the distance from the leading edge of 
the plate to the model, resulted from a joint in the plate. As nearly 
as could be measured, this wave was inclined at the Mach angle) its 
effect on the flow was therefore considered negligible. 

For small values of b/L, a steady flow pattern was obtained 
(figs. 7(a) to 7(c)). For values of b greater than lj inches, how- 
ever, unsteady flow of the type shown in figure 7(d) was obtained. For 
such flows the shock angle could not be measured, but a mean shock angle 
could be calculated from the pressure measurements in the separated 
region. For the steady configurations, the upstream edge of the waves 
that originate near the point of separation is curved in the vicinity of 
the boundary layer but becomes straight farther away from the plate. 

The initial curved portion results from boundary- layer interaction at 
the root of the shock wave, which was ignored in the analysis. In 
analyzing the data, the angle of the straight portion of the shock wave 
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was measured and compared with the shock angle calculated from the 
pressure in the dead-air region. The accuracy of the measured shock 
angle was low, however, because of the thickness of the trace j hence , 
the values of 9 calculated from the pressure measurements were used 
for comparison of experiment with theory . 


Turbulent Boundary Layer 

The analytical and experimental variations of cp with b/l are 
compared in figure 8(a). The value of e was taken equal to 1 for the 
analytical curves. Shown in figure 8(a) are all points obtained with 

the transition strip in place and with the models located 8^- inches from 

■the leading edge of the plate. The experimental points shown in fig- 
ure 8(a) are not directly comparable with the theoretical curves because 
of the effect of reduced span- thickness ratio at large values of b/L. 

In order to provide more nearly comparable data, all values of Cp 
obtained with span-thickness ratios less than 52 were corrected by an 
increment equal to the difference given in figure 5 between the value 
obtained at the actual span- thickness ratio and the value obtained for 
a span-thickness ratio of 32. These corrected values of Cp are compared 
with the theoretical values in figure 8(b). The quantitative agreement 
between the experimental and the theoretical values of cp is quite 
good for values of b/L between 0.050 and 0.125. 

Unsteady flow appeared for values of b/L greater than 0.15, 
which is considerably below the maximum values predicted by theory. 

This unsteadiness may result from the fact that two configurations are 
possible for large bj a modified bow wave could form or the separation 
could extend to the leading edge of the plate. The effective shock 
angles, computed from the measured mean pressure in the dead-air space, 
seem to increase rapidly just before the flow becomes oscillatory, as 
if a detached bow wave were forming, and then to drop to values that 
are lower than some of those obtained with steady flow. The shock angles 
that would be expected if the flow were separating at the leading edge 

of the plate (X = tan -1 b/L) are shown for comparison and the mean pres- 
sures appear to be of the order of magnitude corresponding to such 
separation. However, no definite conclusions can be drawn from these 
data since the nature of the unsteadiness of the flow is unknown, and 
hence the nature of the mean pressure cannot be established. 


For values of b/L less than 0.075 (b/5 less than about 6), 
experimental values of <p remain higher than the theoretical values 
although they, too, should be equal to the Mach angle for b/L =0. 


2 The small angle due to the rate of increase of displacement thickness 
along the plate should be continuous for b » 0 and consequently should 
result in no finite disturbances. 
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The discrepancy between theory and experiment for low values of b/L 
is believed to result from the fact that the effect of the region of 
interaction of the boundary layer with the shock wave, which was 
neglected in the theory, becomes important as b approaches the thick- 
ness of the boundary layer (compare figs. 7(a) and 7(c)). The fact that 
experimental values remain higher than theoretical values for small 
b/L may indicate that the rate of change of momentum thickness in the 
region of shock - boundary-layer interaction is greater than along a 
solid surface. 


Effect of Changes in Shear Profile 

Shown in figure 9 are the values of Cp obtained at Mq = 1.84 when 

the models were moved to 5^ inches from the leading edge with the transi- 

6 1 
tion strip still in place and the values of Cp obtained for L -Q^ inches 

with the transition strip removed. These values are compared with various 
theoretical curves and with the experimental points from figure 8 for L = 

8^ inches with artificial transition. The span- thickness ratios are 
shown for each data point to provide a valid basis for comparison of the 
effects of changes in the initial boundary layer. 

Contrary to expectations, removal of the transition strip increased 
the strength of the shock for each value of b/L. This strength increase 
may indicate that separation originated in a region of transition rather 
than in the laminar region.^ In such transition regions, ©^' is con- 
siderably greater than in either laminar or turbulent regions, so that 
the average shear stress over the dead-air region should be higher, and 
consequently the shock should be stronger, than for completely laminar 

or turbulent flow. When the l/2-inch model was placed 5^ inches from 

the leading edge, on the other hand, the values of Cp for comparable 

span- thickness ratios are lower than for L » 8^ inches. For 

L s 5— inches, the separation region started only about 3 inches from 
4: 

the leading edge of the plate, so that the 3/4-inch region of laminar 
flow ahead of the transition strip could have a more pronounced 

influence on the separation angle than for L = 8— inches. The theo- 
retical curve for transition 3/4- inch from the leading edge, which was 
calculated from the equations in appendix C, agrees very well with the 


3 It seems probable that the separation itself introduces premature trans- 
ition in the upstream laminar layer, so that separation angles corres- 
ponding to laminar flow ahead of the separation point may be observed 
experimentally only for quite low Reynolds "’vribers . 
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experimental point obtained with the largest span- thickness ratio. 

This point, as indicated in figure 5, is probably close to the true 
two-dimensional value. Theoretical curves resulting from the assumption 
of laminar flow on the plate and turbulent shear on the replacement 
wedge (appendix B) were calculated for comparison and indicate that, at 
the Reynolds number of the present tests, the nature of the boundary 
layer ahead of the separation is more significant than the magnitude of 
the assumed rate of increase of momentum thickness along the dead-air 
boundary. In view of the doubt whether separation of the type considered 
herein can occur without inducing transition ahead of the separation, 
these curves are probably only of academic interest. 


DISCUSSION AND EXTENSION 

Flow Patterns for Various Body-Thickness Ratios 

On the basis of the experimental or analytical results presented in 
previous sections, it appears that the flow patterns obtained for a blunt 
body mounted on a flat plate with a supersonic free stream can be divided 
roughly into four types. These types are illustrated in figure 10 for a 
Mach number of approximately 1.9. 

Type 1 (fig. 10(a)). - When the body is of the same order of thick- 
ness as the initial boundary layer, the shock - boundary- layer interaction 
is of primary importance and the analysis of the present paper is 
inapplicable. Data near a Mach number of 1.8 indicate that the body 
should be greater than about 5 times the boundary-layer thickness before 
agreement with the analysis is to be expected. 

Type 2 (fig. 10(b)). - For a certain range of body thicknesses, the^ 
magnitude of which depends greatly on the free-stream Mach number, details 
of the shock - boundary- layer interaction become unimportant, and the 
analysis of this paper is applicable. The assumption that. £ =1 appears 
to yield satisfactory agreement with experimental results in this flow 
region for the Mach number range investigated. For low supersonic Mach 
numbers, this type of flow may not appear at all. 

Type 3 (fig. 10(c)). - As the body thickness approaches the maximum 
value for which the pattern of type 2 can exist, the flow becomes 
unsteady and the shock pattern appears to oscillate between a modified 
bow wave and the straight shock corresponding to separation from the 
leading edge of the plate. The sketch of figure 10(c) shows how the mean 
pressure in the dead-air region can become less than that corresponding 
to separation from the leading edge in this unsteady-flow region. 

Type 4 (fig. 10(d)). - As the body thickness is increased further, 
the angle corresponding to separation at the leading edge approaches the 
maximum for which an attached shock is possible and a configuration 


2183 


me A W 2418 


13 


similar to that shown in figure 10(d) should appear. Except for the 
projection due to the presence of the plate, the shock location in this 
region should be calculable by the method of reference 3. If the plate 
does not extend beyond the foremost point of the detached shock corre- 
sponding to uniform f low, as calculated approximately by the method of 
reference 3, the presence of the plate should have negligible effect on 
the shock form, and the conventional bow wave, normal at the foremost 
point, should reappear. Whether the intermediate stage shown in fig- 
ure 10(d) is stable remains to be established. 


Axially Symmetric Flows 

Discussion . - An extension of the method, whose completion awaits 
more knowledge of shear profiles, is the computation of the effective 
increases in fineness ratio and decreases in drag that can be realized 
by providing an initial boundary layer or wake .ahead of a blunt-nosed 
axially symmetric body. Such a boundary layer or wake can be provided 
by extending a thin rod ahead of the nose (fig. 11) or by suspending a 
small body upstream of the nose to provide a wake. An effective conical 
nose of small weight can in this manner be provided for a blunt body, 
in the event that such a body is desirable for better utilization of 
space. 

The use of dead-air regions to reduce the drag of blunt bodies has 
been suggested before, but no analysis was available to serve as a guide 
in the selection of a method for producing this dead-air region. The 
method shown in figure 11 seems the simplest and is amenable to analysis 
when the required friction data become available. 

For this configuration, the cone angle and consequently the pres- 
sure drag can theoretically be reduced to as small a value as desired 
by increasing the thickness of the boundary layer or wake relative to 
the nose size. If the trend of the data for two-dimensional bodies 
occurs also for axially symmetric bodies, however, the rods required 
may become excesively long before pressure coefficients approaching zero 
can be realized (see fig. 8). Furthermore, as the ratio of the boundary 
layer or wake thickness to the nose radius increases, there is expected 
to be an increase in friction drag almost corresponding to that which 
would be obtained by replacing the projecting rod with separation cone 
by a solid cone. The optimum fineness ratio for the blunt body plus 
separation cone should therefore be almost the same as the optimum 
fineness ratio for a pointed- nose body. Although the length of rod 
required to yield a given fineness ratio cannot as yet be evaluated 
theoretically, the following analysis provides an indication of the 
chief parameters needed to solve the problem in the range of body 
thicknesses for which the flow is steady and the boundary layer is not 
too thick relative to the body. 
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Analysis . - By an analysis for axially symmetric flow similar to 
that used to obtain equation (3) for two-dimensional flow, the equation 
for the case shown in figure 12 is found to be 



1 + h 


0 



Qc 



( 12 ) 


where 


Q 


c 


P u 

c c cos X 
2 c 

Pq u o 


<^>a+5 r 


2it 


pu 

p o u o 





Pb+5, cos X 
1 c 




y dy 
cos X c 


and where subscript c refers to surface values corresponding to the 
separation cone half angle 

On the assumption that © 1 and © are related approximately as 
they would be if a solid cone replaced the separation cone, another 
equation for ©q 2 /© 0 2 can 1 ° e d- erived - and equated to the value of 

equation (12) to determine the cone separation angle. 
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If t is the average shear force per unit area on the cone that 
replaces the dead-air region, 



1 + 


n u 

c 



(13) 


where A c is the surface area of the cone and is given by 


A c = jtb 2 /sin 


f — « 1. From equations (12) and (13), 
"b 


1 + h„ 


U c X 

u^ COS X c 


)■ 


^c 


tan X = T„ tan X 
c c ^ 


o„u_ 0„ 


If the boundary- layer thickness on the projecting rod remains thin 
relative to the rod diameter, then the two-dimensional skin- friction 
formulas should be adequate for calculating 9 q . Thus, for Sq« 


© Q 2 « 2jta 


pu 

P 0 U 0 


0 



k(M 0 )R x ' 1 / 11 cLx 


or 



2itaZ k (Mo) R Z 


-i/n 


Now let p be the ratio of T c to the shear force per unit area on a^ 
flat plate T p and assume, as in the two-dimensional case, that T p is 
approximately equal to the shear stress at x =■ l. Then 

T c = PT p = Pp q u o 2 k(M 0 )R z - 1 / T1 
and equation (14) becomes 

T|— 1 P b b 
^c ' fcan ^c = •q 2 a l 


( 15 ) 
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Because 


l s* L - b ctn ^ 


equation (15) can be written as 



— . 2 : 
p_b 

n 

Tl-l 


T c tan 


2 a _ 

+ p b T c 
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The value of h„ in T„ becomes approximately the two-dimensional 
u ^ 

form factor for 5 q « a since 
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For a given ratio of rod radius to nose radius a/b, equation (16) 
expresses the variation of b/l with cone angle and Mach number when 
the boundary layer on the rod is entirely laminar or turbulent. 
Unfortunately, the magnitude of p is as yet unknown except for the 
case of laminar flow without initial boundary layer. It therefore 
seems advisable to postpone further discussion of equation (16) until 
a theoretical or experimental basis for estimating p becomes available. 


Effect of Nose Shape 

Most of the preceding discussions have assumed that the body has a 
flat nose or at least a sharp corner where the dead-air region terminates. 
For this type of body (wedge or cone) it is to be expected that the 
angle of the nose would have little effect on the shape of the dead-air 
region as long as the wedge or cone angle of the body is considerably 
greater than the wedge or cone angle of the dead-air region. Evidence 
to substantiate this expectation is given in reference 4, where wedges 
with half-angles from 14° to 90° were placed on the floor of a tunnel 
operating at a Mach number of 2.4. It is concluded in reference 4 that 
the foremost shock angle is approximately independent of wedge angle and 
that consequently the separation angle should be chiefly a function of 
Mach number and Reynolds number. 
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For bodies with curved noses ; the point of contact of the dead-air 
boundary with the body should be the point where the slope of the body 
is equal to the slope of the dead-air boundary (see fig. ll) . For the 
particular case of a circular nose, 

b = r cos \ (17) 


so that, for the two-dimensional case, equation (10) becomes 


T tan X 

r cos XT" 



(18) 


which determines the separation angle as a function of the radius of 
the nose. If the nose is noncircular and if the equation of the nose 
contour is 


y = f(x) 


then b is the value of y for which y' = tan 


Supersonic Diffusion by Means of Boundary-Layer Separation 

Because the dead-air region ahead of blunt bodies tends to form a 
wedge or a cone for a range of values of b/L, it is feasible to form 
a supersonic inlet by placing a cowling around a configuration such as 
that shown in figure 11. Such an inlet should yield pressure recoveries 
that differ little from those obtained with inlets using solid cones or 
wedges as long as the mass flow into the inlet is the maximum possible. 
When the mass flow is reduced, this type of inlet may have lower 
additive drag due to spillage than the solid-body inlet, because the 
dead-air region will not support the pressure gradient across a strong 
detached shock wave. Consequently, the size and angle of the dead-air 
region should change as the mass flow is reduced. It remains to be 
demonstrated whether the separation inlet can be stable under reduced 
mass flow conditions. 


CONCLUDING REMARKS 


An analysis of supersonic flow against blunt bodies located in 
boundary layers or wakes indicates that wedge-type or conical-type dead- 
air regions occur only over a limited range of body thicknesses relative 
to the initial boundary layer or wake thickness. For the particular case 
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of two-dimensional Hunt bodies mounted on a flat plate, quantitative 
agreement between predicted and experimental wedge- separation angles 
was obtained in this range of body thicknesses when the change in 
momentum thickness along the dead-air region was assumed to be equal to 
that which would be. obtained if a solid wedge replaced the dead-air 
region. When the body thickness approached the maximum for which wedge- 
type separation was possible, unsteady flow was observed. For body 
thicknesses of the same order as the initial boundary layer thickness, 
shock angles larger than those predicted were obtained. This discrep- 
ancy was attributed to neglect of boundary- layer - shock-wave interaction 
in the analysis. 

For low supersonic Mach numbers, the two-dimensional analysis 
indicated that the simple wedge-shaped dead-air region can exist only 
for body thicknesses that are of the same order of thickness as the 
initial boundary layer. For such bodies, however, details of the shock - 
boundary- layer interaction, which were neglected in the analysis, became 
important. The analysis given herein should therefore be of interest 
chiefly at high Mach numbers, where the cone or wedge-type separation 
can exist for bodies that are relatively thick compared with the initial 
boundary layer. 

The unsteady flow observed when the body thickness approached the 
maximum possible thickness for wedge-type separation may result because 
two types of configuration are possible. A detached bow shock could 
form or the separation could extend to the origin of the boundary layer. 
Experimental values of mean pressures near the nose of the body indicate 
that the latter configuration may predominate for body thicknesses 
sufficiently large relative to the length of the plate upstream of the 
body, but no steady configurations of either type were observed. 

The analysis developed for wedge-type dead-air regions was also 
applied to conical dead-air regions ahead of axially symmetric bodies, 
but for this case quantitative estimates must await theoretical or 
experimental determination of the skin friction on cones in the pressure 
of an initial boundary layer at supersonic speeds. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, May 2, 1951. 
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APPENDIX A 


SYMBOLS 

The following symbols are used in this report: 


area 


radius of rod ahead of axially symmetric blunt body 

thickness of two-dimensional body or radius of axially 
symmetric body at point of contact with dead-air region 

form factor , S */© 

-l/n . , 

function of Mq which when multiplied by R yields 

friction drag coefficient 

distance in x-direction from start of boundary layer to point 
of body where thickness or radius is b 

distance from start of boundary layer to beginning of dead-air 
region 

Mach number 


power-law profile parameter for turbulent boundary layer, 
assumed equal to 7 


pressure coefficient, 


Prandtl number 



static pressure 

function of Mq and X defined by equation (3) 
function of Mq and \ c defined by equation (12) 
Reynolds number, based on distance given by subscripts 
radius of nose of circular-nosed bodies 
distance normal to cone or wedge surface 
function of Mq and X defined by equation (6) 
function of Mq and defined by equation (14) 
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u 

local velocity 


X 

distance in free- stream direction 


y 

distance normal to free- stream direction 


r 

ratio of specific heats 


5 

boundary- layer thickness 


5* 

displacement thickness of boundary layer 


e 

factor of proportionality between change in e on flat 
plate and along dead-air region 


n 

power parameter in skin- friction laws 


9 

shock angle resulting from dead-air region 


\ 

angle of boundary of dead-air region 

• 

T 

shear stress 


9 

momentum thickness of boundary layer 


P 

density 


P 

ratio of skin friction on cone to skin friction on wedge 


Subscripts: 


av 

average 


0 

conditions ahead of dead-air region 


1 

conditions downstream of shock or at downstream end of 
dead-air region 


c 

conditions on surface of cone used to replace conical 
dead-air region 


P 

flat -plate value 


Za 

laminar value 


max 

maximum 


t 

turbulent value 



Primes denote differentiation with respect to x 
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APPENDIX B 


LAMINAR BOUNDARY LAYER WITH TURBULENT SHEAR ON 
SEPARATION WEDGE 

The general equation for evaluating the separation angle is 
( see equation ( 6) ) : 


0 ' 

€ —— b s« T tan A. (Bl) 

9 0 


It is desired to compare the values of b obtained for a given T 
and \ (or Mq and X) under two assumptions: 

(l) Laminar flow on plate and on replacing wedge 


(2) Laminar flow on plate and turbulent flow on replacing wedge 


If values corresponding to assumptions (l) and (2) are identified with 
subscripts 1 and 2, respectively, the following relation applies for 
given Mq and \ : 



b 2 (V^O^ 
= ( e oW 2 


(B2) 


Now (e A '/e ri ) = — — from equation (8) of the text. In the quantity 
u 

(9q'/©q) , 0q is the momentum thickness corresponding to a laminar run 
2 

of distance whereas 0 q' is the turbulent friction coefficient 

corresponding to the momentum thickness © . Thus, 




2Z 2 k 7a R 2 


- 1/2 

2 



-1/7 
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where Z 2 ' is the equivalent turbulent run required to build up the 
actual momentum thickness From these equations } 


and 


or 



V 


12 

7 




k t R 2, 


- 1/2 

^7? 53 



l 


2 


12 k 7a _ -5/14 

J ^ h 


7/6 


Equation (B2) now becomes 

yr - (v) 1/6 (v V 5/14 ) ?/6 ■ \ 5 , 1 T <B3) 

To estimate the magnitude of the difference resulting from cases (l) and 
(2), assume that k^/k^. is nearly independent of Mach number and has 

the incompressible value^ 



1.528 

0.0131 


101.3 « 100 


Then 


b 2/ Z 2 _ 255 

h 7 h ' Rj 2 °- 416 


(B4) 


^The value of k^aAt obtained for Mq * 1.9 from equations in refer- 
ences 1 and 2 is 103.2 for the Reynolds number of the tests reported 
herein. 
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If the models are placed at the distance L from the leading edge of 
the plate, equation (B4) becomes 


b 2 /L 




-,0.584 


\ 


235 


0.416 


R, 


-r=- ctn \ 
Jj 


(B5) 


From equation (B5) , b^/L can "be calculated from b^/L for any value 
of M q and X. As might he expected, the variation of b 2 /L with X 

or cp is no longer independent of Reynolds number. Values of b/l 
obtained from equation (B5) for two values of L are compared with the 
purely laminar values in figure 9. 
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APPENDIX C 


COMPUTATION OF SEPARATION ANGLE WHEN UPSTREAM BOUNDARY 
LAYER IS PARTLY LAMINAR 

The extent of the laminar region "before transition influences the 
separation angle through the magnitude of 0 Q , which is proportional o 
the integrated friction drag up to the separation point. The separation 
angle can he calculated as function of b/L for this case if the 
transition is assumed to occur at a definite point; x - an 1 
value of 9 at this point can be estimated. For the present purpose. 

It is assumed that the jump in 0 at l ± is one-half the difference 

between the value of 0 obtained with the turbulent equation and the 
value of 0 obtained with the laminar equation. This assumption 
represents a compromise between two conventional assumptions, namely, 
that there is no discontinuity in © at a transition point or that the 

“ 9 at the transition point is such that the values of 9 beyond 

the transition point are the same as they would be with turbulent flow 
from the leading edge. The relative magnitudes of 0 Q resulting from 
these two assumptions are indicated in the following sketch. The com- 
promise values used should be adequate for estimating the magnitude of 
the effect of a partly laminar boundary layer. 


Equation 



With this assumption the equation for 0 Q can be written: 


n- 


P 


9 0 " 


- 1/2 


k la R x ' + \ h KV 1/7 ■ k Jh R J 1 ) * 


k tV 1/7 


- »1 (| k laV l/2 - I k t R l 


1 1/? ) + i WH 


-1/7 
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Since 

V - *t E z 2 ' 1/7 


equation (6) in the text becomes 


9'h 

=s T tan \ ss 

9^ 


cb *ty 1/7 


\ I 2 Vli 


-i/2 _ 2 


- k, R, _1 / 7 \ + 1 7_k. R, 


3 


6 2 t J, 


-1/7 


(Cl) 


To compare the values of h obtained from equation (Cl) with those 
obtained for completely turbulent flow, let T tan X (or Mq and X ) 
remain the same for both cases. If values corresponding to completely 
turbulent flow are denoted by the subscript a and those for partly 
laminar flow by subscript (3, 


(C2) 



After some algebraic manipulation, equations (Cl) and (C2) result in the 
following equation: 


, \ 6 / 7 


(b/Z2)0 


'iso Rjj 5 / 14 - §y \ 



!(t 


(C3) 


a 


where k^/k^ has been taken equal to 100, as in appendix B. Because 
2 2 =5 L - b ctn equation (3) can also be written 

(t>/L) E E ~(l! CtD I 7 ' (b/L)„, 


/ 1 6/7 r 


6 / 7 , ( 

b\ 

(t) ( 150 e Ii ' 5/14 ' l) + ' 

• (4 ctn 

X -{ 

L A 


ctn X 
(C4) 

This equation can be solved for (b/L)^ by trial when Z-, and the 
Reynolds number are known since (b/L) a as a function of \ is known 
from computations for the case of completely turbulent boundary layers. 
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In figure 9, the variation of b/L with CP from equation (C4) for a 
3^4-inch laminar run before transition is compared with the variation 
for completely turbulent flow and with experimental results. 
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(a) Supersonic initial velocities. (t>) Partly subsonic initial velocity. 

Figure 1. - Flow against blunt bodies in nonunifomi stagnation pressure field. 



Figure 2. - Sketch used for analysis of flow against blunt bodies mounted 

on flat plate. 


Form factor 



Figure 3. - Form factor for laminar and turbulent compressible, flat-plate boundary- 
layer profiles as function of free-stream Mach number. 
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Separation angle, A, , deg 
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Body thickness ratio, h/L 

(a) Separation angle. 

Figure 4. - Variation of separation angle and shock angle with body thickness ratio. Two-dimensional flow. 
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Turbulent (t) = N = 7) 

— T^mlriRT 



Body thickness ratio, b/L 
(b) Shock angle. 

Figure 4. - Concluded. Variation of separation angle and shock angle with body thickness ratio. Two-dimensional flow, 
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Figure 5. - Effect of span-thickness ratio on separation parameters. 
Free -stream Mach number Mq, 1.84; b, 0.5 inch; L, 5.25 inches. 
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(b) Ratio of 6 0 ' to 9 0 . 


Figure 

used 


6 - Experimental and empirical boundary -layer parameters for flat plate 

in tests. Artificial transition; free-stream Mach number M Q , l.bB. 
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(a) Free-stream Mach number Mq, 1.84; 
b, -i inch; L, 8rj inches; artificial 
trans ition. 



(c) Free-stream Mach number Mq, 1.84; 
b, 1 inch; L, inches; artificial 
transition. 


(b) Free-stream Mach number Mq, 1.84; 
b, inch; L, 8^ inches; artificial 
transition. 



(d) Free-stream Mach number Mq, 1.84; 
b, 1-jg inches; L, 8^ inches; artificial 
transition. 



1 


(e) Free-stream Mach number Mq, 1.84; 
b, -i inch; L, 5^ inches; artificial 
transition. 


( f ) Free-stream Mach number Mq, 1.87; 
b, inch; L, 8^ inches; natural 
transition. 

C- 27618 



Figure 7. - Schlieren photographs of flow against rectangular rods mounted on flat plate. 




NACA TN 2418 


35 



(a) Uncorrected data. 


Figure 8. - Comparison of experimental and theoretical variation of separation shock angle with 
body thickness ratio. L, 8^ inches; artificial transition. 



Shock angle, 
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(b) Corrected data. 


Figure 8. - Continued. Comparison of experimental and theoretical variation of separation 
shock angle with body thickness ratio. L, 8i inches; artificial transition. 
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Body thickness ratio, b/L 


Figure 9. - Effect of changes in boundary layer ahead of separation point. 
Free -stream Mach number Mq, 1.84. 
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Figure 10. - Observed and probable flow regimes for various ranges of body 
thickness ratio. Free -stream Mach number Mq, approximately 1.9. 
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Figure 11. - Reduction of drag of thick bodies by provision of initial boundary 

layer . 



Figure 12. - Sketch for analysis of axially symmetric flow against blunt 
bodies with initial boundary layer. 
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